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An introduction to tiie tiieory of critical behavior at Lifshitz points is given, and the recent 
progress made in applying the field-theoretic renormalization group (RG) approach to (j)'^ n- 
vector models representing universality classes of m-axial Lifshitz points is surveyed. The 
origins of the difficulties that had hindered a full two-loop RG analysis near the upper critical 
dimension for more than 20 years and produced long-standing contradictory e-expansion re- 
sults are discussed. It is outlined how to cope with them. The pivotal role the considered class 
of continuum models might play in a systematic investigation of anisotropic scale invariance 
within the context of thermal equilibrium systems is emphasized. This could shed light on the 
question of whether anisotropic scale invariance implies an even larger invariance, as recently 
claimed in the literature. 
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1 Preliminary remarks 

The aim of this article is to give a brief introduction to the field of critical behavior at Lifshitz 
points [1,2], and to survey the recent progress that has been made in applying the field-theoretic 
renormalization group (RG) approach [3, 4, 5, 6, 7, 8]. As is appropriate for a conference on the 
topic 'renormalization group', 1 shall assume some basic familiarity of the reader with RG ideas. 
However, in view of the mixed background of the participants of the conference, no extensive 
knowledge of the relevant condensed matter physics is presupposed. 

Since the literature on Lifshitz points — or, more generally, on systems with spatially mod- 
ulated phases — is vast, it is impossible to mention or even cite all relevant papers. The paper 
is designed to give a reasonably self-contained account of the issues on which we focus, and 
to serve as a guide to the literature. The choice of the references has been made accordingly. 
There exist extensive review articles [1] and [2], which survey the literature till 1992 and contain 
extensive lists of references. The reader may consult these for further information on topics that 
had to be left out here. 
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2 Introduction and background 



2.1 Generic phase diagram with a Lifshitz point 



The concept of a Lifshitz point was introduced more than 25 years ago [9]. It is a point in the 
phase diagram at which a disordered phase, a spatially homogeneous ordered phase, and a spa- 
tially modular ordered phase meet. A typical phase diagram with a Lifshitz point is depicted in 
Fig. L In the case of a ferromagnet, the disordered and uniform ordered phases are the usual 
paramagnetic and ferromagnetic phases. The order of the latter corresponds to an infinite 'mod- 
ulation' wave-length, i.e., a modulation wave-vector Qq = 0. The transition that occurs upon 
crossing the phase boundary between the disordered and the uniform ordered phases is contin- 
uous ('second-order transition'). Thus each point on the line T = Tc{g) separating these two 
phases and emerging from the Lifshitz point is a critical point. The variable <? is a second (in- 
tensive) thermodynamic variable besides temperature T. What it stands for depends on the type 
of system considered: In the case of organic crystals like TTF-TCNQ [10, 11], g corresponds 
to pressure; in the cases of the much studied magnet MnP [12, 13] and the so-called ANNNI 
model [14], which we will both briefly consider below, g stands for the magnetic field compo- 
nent perpendicular to the order parameter and a ratio of an antiferromagnetic to a ferromagnetic 
interaction coefficient, respectively. The important point to remember is that g does not couple 
directly to the order parameter. (If it did, a small change g ^ g -\- 5g,T ^ T^-\- along any 
direction in the gT plane — and specifically along Tc(.9) — would destroy the critical behavior; 
i.e., there would not be a critical line.) For this reason, g is commonly called a nonordering 
thermodynamic field. 



Fig. 1 . Schematic phase diagram with a Lifshitz point L. The disordered phase is separated from the uniform 
and modulated ordered phases via a critical line. The crossover exponent ip is defined via the behavior of 
the critical line near L; the wave- vector exponent fiq describes how the modulation wave-vector go tends 
to the value go,L as L is approached along the critical line between the disordered and modulated ordered 
phases. 
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The term 'thermodynamic field' is probably not too familiar among high-energy physicist, 
and should by no means be confused with what is meant by a field in field theory. Its usage 
was suggested in a seminal paper [15] by Griffiths and Wheeler for intensive thermodynamic 
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variables like T, magnetic field, pressure, and chemical potential that take the same values on 
both sides of ( in general, first-order) bulk phase transition. Their thermodynamic counterparts 
are the densities of extensive variables like the magnetization, energy, and entropy densities that 
normally jump at the transition (though not necessarily all of them) since their values in the 
respective pure bulk phases differ. 

In the modulated ordered phase, the order is characterized by a nonzero modulation wave 
vector Qq, which varies with T and q. The transition between the disordered and modulated 
ordered phase is continuous as well. Hence the Lifshitz point divides the critical hne Tc{g) into 
two sections. The transition between the homogeneous and modulated ordered phases can be of 
first or second order; for models with a scalar order parameter it is generically discontinuous, for 
specific models with vector order parameters it is found to be continuous [16]. 

Let us also note that the modulation wave- vector does not necessarily have to vanish at the 
Lifshitz point. To see this, recall that in the case of antiferromagnets, the order parameter which 
acquires a nonzero value in the uniform ordered phase is not the magnetization but the staggered 
magnetization. If we consider lattice models, this is given by the average of the sum of all spins 
on one (up-spin) sublattice minus the sum of those on the other (down-spin) sublattice. Thus the 
homogeneous ordered phase does not correspond to a homogeneous magnetization density but 
to one modulated with the corresponding nonzero wave- vector at the boundary of the Brillouin 
zone. This value of the modulation wave-vector <7q applies in particular to the Lifshitz point, so 
9o L 7^ 0. For the sake of simplicity, we shall always use ferromagnetic language in the sequel, 
taking Qq.l to vanish. 

2.2 Critical exponents and continuum models 

To reach the Lifshitz point, both the temperature T and the nonordering field g must be fine- 
tuned. This tells us that g must correspond to a relevant variable in the RG sense. We denote 
the associated crossover exponent as ip. It describes the behavior of the critical line T^i^g) in the 
vicinity of L: As indicated in Fig. 1, we have 

5T, = T,{g)-nr^\5g\^l'P , Sg = g - gi^ . (1) 

In order to describe the behavior of the modulation wave-vector near L, one introduces a wave- 
vector exponent (3q via 

6qo = qo- go,L ~ {Sgf- . (2) 

There are further critical exponents that are needed to characterize the critical behavior at Lifshitz 
points. We will introduce some of these below. 

How a Lifshitz point can occur can be easily understood within Landau theory. Landau 
theory leads one to consider the following natural generalization of the usual (j)'^ model with the 
Hamiltonian (Euclidean action) 

n = Jd'^x[^iv^<pf + f (V||0)^ + f (A||0)^ + "i^' + ^\'i>\']- (3) 

Here (f) — (0") is an n-component order-parameter field. The d-dimensional position vec- 
tor X = (a;|| , has an m-dimensional 'paraUel' component a;|| and a {d — m)-dimensional 
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'perpendicular' one x±. The coefficients of the squared gradient terms, and those of the Hamil- 
tonian's other interaction terms, generally depend on the thermodynamic variables T and g. We 
have assumed that the squared gradient terms involve just two distinct coefficients, where the 
one of (V^0)^ remains positive and has been transformed to unity by means of an appropriate 
choice of the amphtude of <p. Its parallel counterpart, po, is permitted to change sign. To ensure 
stability when po < 0, the term (A||(/>)^ with a positive coefficient ao > has been added. 
The critical line between the disordered and homogeneously ordered phases is given in Landau 
theory by tq = with po > 0. The Lifshitz point is located at tq = po = in this classical 
approximation. 

Near the Lifshitz point the interaction coefficients can be expanded about T = Tl and g = 
5l- For the coefficients that remain positive, i.e. uq and ctq, the expansions may be truncated at 
zeroth order, so that uq and ctq become independent of T and g. The deviations Stq = tq — ro,L 
and Sao = o-q ~ ctq.l change sign at L (where tq.l and (To,l vanish in Landau theory); their 
expansions must be retained to linear order in T — Tl and 6g. 

3 What is interesting about studying critical behavior at Lifshitz points 

One reason for the ongoing interest in critical behavior at Lifshitz points is the wealth of distinct 
physical systems with such multi-critical points. They range from magnetic ones [12, 13, 17], 
ferroelectric crystals [18], and charge-transfer salts [10, 11] to liquid crystals [19], systems un- 
dergoing structural phase transitions [20] or having domain-wall instabiUties [21], and the so- 
called ANNNl model [14]. Lifshitz points have been discusssed recently even in the context of 
superconductors [22] and polymer blends [23]. 

From a more general perspective, the problem is interesting because it provides well-defined 
clear examples of systems exhibiting anisotropic scale invariance (ASI). At conventional criti- 
cal points, scaling operators O such as the order parameter and the energy density transform 
asymptotically as 0{£x) = 0{x) under scale transformations, with a (in general nontriv- 
ial) scaling dimension Aq- In the case of ASI, the position coordinates x, or the position and 
time coordinates (in the case of time-dependent phenomena), divide into two (or more) groups, 
say X = (a;||, that must be scaled with different powers of the scale factor i to recover 
0{x); one has 

0{e^xii , ex_L) = 0{x) , (4) 

with an anisotropy exponent 6 different from unity. That ASI appUes at Lifshitz points is evident 
from the Hamiltonian (3): The momentum-space two-point vertex function of the free theory 
obviously behaves at the Lifshitz point tq = po = as 

f(2)(q||= 0,9^^0)^51-"- and f (q,, ^ 0, = 0) ^ gp" , (5) 

where the analogs 77L2 and 77L4 of the usual correlation exponent 77 take their mean-field value 
zero, but differ from it beyond Landau theory. In terms of these, the anisotropy exponent reads 
6* = (2 — ?7l2)/(4 — 1JL4:). The asymptotic validity of ASI means that two distinct correlation 
lengths and ^[| are needed to characterize the region within which (f>{x) is correlated to 4>(0); 
these diverge at 5 = 5l as functions of ST = {T - Tl)/Tl like {STl-"^^ and |(5T|-'^" with 
different exponents 1/1,2 and !/l4 = 9 ui,2, respectively. 



Critical behavior at m-axial Lifshitz points 



5 



A familiar arena of ASI are dynamic critical phenomena near equilibrium [24, 25]. In 
their case the time t and position x must be scaled differently; the analog of Eq. (4) becomes 
^{tx, i^t) = t~^'i' 4>{x, t) for the order parameter, and the role of 9 is played by the dynamic ex- 
ponent 3. Their simplifying feature is that detailed balance and fluctuation-dissipation theorems 
hold. This entails that the stationary states of the dynamics are guaranteed to be thermal equilib- 
rium ones described by an a priori known Hamiltonian. Hence the problem of the steady-state 
correlations splits off from the dynamics; all critical exponents of static origin can be determined 
without dealing with dynamics, only genuine dynamic properties require the analysis of the dy- 
namic field theory. 

ASI is abundant in non-equilibrium systems such as driven diffusive systems [26] and surface 
growth processes [27], which have attracted considerable attention during the past decade. In 
their case detailed balance is not normally valid. Finding the steady-state solutions therefore is a 
nontrivial task, requiring the investigation of the long-time limit of the dynamics. This must be 
solved before the ASI of the steady-state correlations can be investigated. 

Two features make critical behavior at m-axial Lifshitz points a very attractive stage for 
the study of ASI: (i) One can stay entirely within the realm of equilibrium statistical physics, 
and (ii) the class of models (3) involves a parameter, m, that can be varied. Some time ago 
Henkel [28, 29] argued that ASI should imply additional invariances, just as scale and rotational 
invariance in local field theories lead to the larger symmetry group of conformal transformations 
[30, 31, 32]. He made predictions for the form of the scaling function of the pair correlation 
function at criticality, which are in conformity with analytic results for certain spherical models 
[33]. Recent Monte Carlo results for an equilibrium [34] and non-equilibrium systems [35] 
appear to support these claims. Yet the general validity of such additional invariances has not 
been shown, nor is there a good understanding of their origin or the conditions under which they 
hold. For example, in Ref. [28] the assumption that the anisotropy exponent takes the rational 
values 2/p, p e N, is needed to ensure that the considered sub-algebra closes. However, in 
the case of the uniaxial Lifshitz point in d = 3 dimensions studied by means of simulations [34], 
this condition is unlikely to be fulfilled since the e expansion [5, 6] reveals that 9 differs from 
1/2 at order e^, although the d = ?> estimate 9 ^ 0.487 it yields for m = n = 1 is pretty close 
to it. Building on the field-theoretic analysis described in Refs. [5] and [6], and outlined below, 
one should be able to clarify whether the suggested 'generalized conformal invariance' and the 
predicted form of two-point scaling functions hold indeed. 

Before turning to the RG analysis of the models (3), let us discuss two examples of systems 
with a Lifshitz point. 



4 Two examples of systems with Lifshitz points: the ANNNI model and MnP 

We begin with the axial-next-nearest-neighbor Ising (ANNNI) model, defined through the Hamil- 
tonian 

HaNNNI = - Si Si+s^ - ^ '^i+h - Y ' 

i,(5|| ?,i5|| 

where Sj ± 1 are Ising spins residing on the sites i of the cubic lattice illustrated in Fig. 2. Here 
5\\ and 8',, are nearest-neighbor (nn) and next-nearest-neighbor (mm) displacements along one 
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axis, while 5± denote nn displacements along any of the other (perpendicular) directions. The 
nn couplings are ferromagnetic, Jq > 0, Ji > 0; the nnn bond is antiferromagnetic, J2 < 0. 
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Fig. 2. ANNNI model. Along one axis the spins are coupled via ferromagnetic nearest-neighbor bonds of 
strength Ji and antiferromagnetic axial next-nearest neighbor bonds of strength IJ2I = — J2. Along the 
other directions there are only ferromagnetic nearest-neighbor interactions of strength Jo. 

The phase diagram of the three-dimensional ANNNI model is depicted in Fig. 3. Obviously, 
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Fig. 3. Phase diagram of the three-dimensional ANNNI model, taken from Ref. [2] (with permission kindly 
granted by its author); to illustrate the ordering in the (2,2) phase the inset with the arrows has been added. 
L is a uniaxial Lifshitz point. 

the nonordering field g here translates into the ratio k = —J2/Ji- The phase labeled (2,2), also 
called (2) structure, corresponds to an antiferromagnetic ordering of the kind indicated by the 
arrows in the inset. The region marked 'modulated' actually has considerably more structure: 
From the multi-phase point fceT/ Ji — 0, k = 1/2, commensurate phases of type (3) and (2^3), 
p — 1,2,... , 00, split off. Here (3) denotes a periodic layer sequence of 3-bands, while (2^3) 
signifies a periodic layer sequence of p 2-bands followed by a 3-band, where a fc-band means 
that the magnetization has the same sign in k successive layers. Though interesting, these details 
cannot be expounded here because they are off our main topic. They can be found in the literature 
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[2, 14]. The essential point for us to note is that L is a M«/ax/flZ (m = 1) Lifshitz point. It has been 
repeatedly studied via Monte Carlo simulations to determine the values of the critical exponents 
for the case m = n=l,d = 3 [36, 37, 38, 34]. 

An experimentally much studied system is the orthorhombic metallic compound MnP [12, 
13, 39]. Its phase diagram differs, depending on whether the magnetic field H is directed along 
the a, b, or c crystal axis. Both for H\\a and H\\b a uniaxial Lifshitz point is found, but not for 
-Ffjlc. The phase diagram of Ref. [12] for the case H\\bis reproduced in Fig. 4 (with permission 
kindly granted by the first author, Y. Shapira ). The magnetic field component H ■ b plays the 




Fig. 4. Phase diagram of MnP for a magnetic field parallel to the b axis, according to Refs. [12] and [13]. 

role of the thermodynamic field g. In the phase labeled 'para' the magnetic moments ('spins') 
s{x) are aligned along the H oc b direction. In the 'ferro' phase the magnetization m = {s) has 
a component along the c axis and hence is tilted with respect to b. The component c • s plays the 
role of the order parameter; since a is a very hard axis, s ■ a can be ignored [39]. In the fan phase 
the spins rotate in the be plane (but make no full turn as in the screw phase 'SCR' not considered 
here); there is modulated order with a modulation wave-vector Qq oc a. The ferro-fan transition 
is first order. The meeting point of the para, ferro, and fan phases is an m = n = 1 Lifshitz 
point. 



5 Dimensionality expansions and renormalization group analysis 

Setting po = To = in the Hamiltonian (3), we see that q± scales as g| at the Lifshitz point of 
the Gaussian {uo=0) theory. By dimensional analysis we have [x±] = [ary] = a-y^/i"^/^, 
[to] = /U^, [po] = ctq^^/U, and [uo] = <''o*^^M^ with e = d*{m) — d, where 

d*{m)=4+—, m<8. (7) 

From the given engineering dimensions we can read off how the interaction coefficients transform 
under the scale transformation — > ij,£. The coupling constant uq becomes marginal at d = 
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d*{m), the upper critical dimension (UCD). Note that m = 8 implies rf* = m = 8; i.e., this 
is the case of the isotropic Lifshitz point in which only the parallel components of position and 
momentum remain. 

One can analytically continue the parallel and perpendicular momentum integrals in their 
respective dimensions m and d — m, and hence take those as continuously varying. Ideally 
one would like to determine the expansions in m and d — m of the critical exponents and other 
universal quantities about any point of the line d*(m). The situation is illustrated in Fig. 5. The 




1 2 3 4 5 6 7 8 



Fig. 5. md plane with the line of upper critical dimensions d* (m). For further explanations see main text. 

point (m, d) = (0, 4) marks the UCD of the standard (isotropic) l^l"* model; analyzing it via the 
conventional e — 4 — d expansion means to move away from this point along the path indicated 
by J|. More generally, e expansions at fixed m correspond to paths parallel to the d-axis. How 
to reach the physically interesting point (1, 3) along such a path is indicated. In the case of the 
isotropic (m=rf) Lifshitz point, things are different: The conventional expansion in eg = 8 — 
is along the diagonal to = d [8, 9]. 

We must also keep in mind an obvious condition points (m, d) to which extrapolations make 
sense must satisfy: d must be larger than d^:{m, n), the lower critical dimension (LCD) below 
which a Lifshitz — or, if m = 0, a critical — point cannot occur. Since the LCD of the Ising model 
with short-range interactions is d*(0, 1) = 1, we clearly must have d > 1. In the n > 1 case 
with the continuous 0{n) symmetry, familiar spin-wave arguments and an analysis based on an 
adequate nonlinear sigma model [40] show that the ordered phase becomes thermally unstable at 
any T > if < d?^"^ (to) = 2 + f (cf. Fig. 5). 

The existence of a Lifshitz point with < to < 8 also requires that the modulated ordered 
phase be thermally stable. In the critical region of this phase, fluctuations of the Fourier com- 
ponents (f)q of the order parameter with wave- vectors q ~ ±qg are dominant. As pointed out 
in Ref. [41], one therefore expects that an n-component system with a heUcal structure behaves 
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critically as a 2n-component (f)"^ model whose cf)'^ terms are 0(ri) but not 0{2n) symmetric. 
Specifically for n = 1, one arrives at an anisotropic two-component model. It has (besides 
others) an 0(2) symmetric fixed point that is believed to be stable [41]. This suggests that the 
long-range order (LRO) of the modulated phase should be destroyed in d < 2 dimensions by 
thermal fluctuations at any temperature T > 0. Furthermore, if the presumed isotropy of the 
Hamiltonian (3) in the parallel subspace can indeed be taken for granted, one can exploit the in- 
variance under arbitrary rotations along the lines of Mermin and Wagner [42] to show the absence 
of a helical phase with orientational LRO form < d < m + 1 [43, 44].^ 

The goal of expanding in m and d about a general point on the line d*{m) was envisaged 
already in Ref. [9]. Yet its realization turned out to be extremely difficult. Computing the e ex- 
pansions of the critical exponents for general values of m to first order in e is easy. Their 0(e) 
coefficients are independent of m, and essentially determined by combinatorial factors. (The op- 
erator product expansion is structurally similar to the one that applies to the standard 1(6 1** theory. 
Upon generalizing Cardy's analysis in Ref. [45] to the m > case, one can determine these 
coefficients without having to work out the Feynman diagrams in detail.) The real challenges 
start at order e^. 

The first results to order for general m were given in 1 977 by Mukamel [46] . The term of 
i%2 he determined by means of Wilson's momentum- shell integration method was independent 
of m, that of ?7l4 differed from it by a simple A/i-dependent factor. He also computed /3q to 
O(e^) for m < 6. For m = 1, these results were confirmed by Hornreich and Bruce [47]. 
Utilizing also Wilson's technique, Sak and Grest [48] did an independent calculation; because 
of the severe technical difficulties they encountered for general m they confined themselves to 
m = 2 and m = 6. Their results are at variance with Mukamel's; in particular, the coefficients 
of 77L2 (w, n) for m = 2 and m = 6 are not equal, and hence not m-independent. 

The apphcation of modern field-theory RG approaches to the problem began in 1998. Mer- 
gulhao and Carneiro formulated normalization conditions and derived RG equations for the 
renormalized theory [3]. In a subsequent paper [4], they reproduced Sak and Grest's results 
for 77l2 and ?7l4 with m = 2 and m = 6, and performed a two-loop calculation for these two 
values of m. They fixed the perpendicular dimension d — m at d* (2) — 2 = 3 and d* (6) — 6 = 1, 
taking the parallel one as 2 — ey and 6 — e||, respectively. The respective paths starting from the 
points (m, d) = (2, 5) and (6, 7) are indicated in Fig. 5. 

In two recent papers [5, 6] Shpot and myself have been able to perform a full two-loop 
calculation for general m € (0, 8) and to determine the e expansions of all critical, crossover, 
wave- vector, and correction-to-scaling (Wegner) exponents to O(e^). The results are analytical 
except that the two-loop terms of the required renormalization functions and the series expansion 
coefficients of the exponents' terms involve four well-defined single integrals j^{m), jaimn), 
jp{m), and Ju{fn) which for general m we have not been able to evaluate analytically, though 
for the special values m = 0,2, 6, and 8. For other values of m theses integrals can be computed 
numerically, as we did for to = 1, 2, . . . , 7. 

These results stand a number of nontrivial checks. First of all, they reduce to the well-known 



^In Ref. [44] the fact that the one-loop shift of Tc{g) on the helicoidal section of the critical line diverges in the 
infrared \fm < d < m + 1 is interpreted as signaling the absence of helical LRO. However, this property alone is not 
sufficient to rule out LRO: The corresponding shift of of the standard one-component model is infrared divergent 
for d<2. Nevertheless, the d = 2 Ising model has a ferromagnetic low-temperature phase. 
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results for the standard \<p\'^ model in the limit m 0.'^ Second, the analytical results they 
yield for m = 2 and to = 6 are consistent with and extend those of Sak and Grest [48] and 
of MerguUiao and Cameiro [4]. [The original O(e^) results of Ref. [4] for vi,2im=2, 6) and 
t'L4(™=2, 6) disagreed with ours but become identical to those upon elimination of two minor 
computational errors.] The isotropic case m = d provides a third, highly nontrivial, check. To 
see this, note that the e expansions of the critical exponents A = r}u, ... is of the form 

A(n,TO,d)|^^^.(„)_^ = Ao + Ai(n)e + A2(n,TO)e2 + 0(e3) , (8) 

i.e., the m-dependence starts at order e^. To compare with the es = 8 — d expansion for the 
isotropic case, we set m = d = 8— eg (which gives e = e8/2)inEq. (8). Hence the expansions of 
the critical exponents A of the isotropic Lifshitz point to quadratic order in es should result from 
the right-hand side via the replacements e es,/2 and X^in, to) A2(n, 8—). The so-obtained 
eg expansions of the exponents A = ?7L4(^^, d, d), 1^4 (n, rf, d), (p{n, d, d), (3q{n, d, d), and the 
correction-to-scaling exponent LUi^^in, rf, d) have been verified by means of a direct field-theo- 
retic investigation of the to = d case [8]. Yet, there is one group of authors [49] who obtained — 
and still favor — results at variance with ours in Refs. [5, 6]; their findings and criticism have been 
refuted in Refs. [7, 8]. 

The origin of the technical difficulties which had prevented two-loop calculations for so long 
and caused the mentioned long-standing controversies can be traced back to the difficult form of 

the free propagator G{x) at the Lifshitz point tq = po = 0. In the isotropic cases to, = and 
m = d, G{x) is a simple power of x. However, for general m (and cto = 1) it becomes 

G{x\\,x±)= / 4 I 2 = xl^+^^m,d{x\\^± ) ' (9) 

J q 111+ 11 

where ^rn,d are extremely complicated scaling functions of the form 

The latter are generalizations of generalized hypergeometric functions known as Fox-Wright i^'i 
functions [5, 6]. They have asymptotic expansions in powers of v~^. On the line d*(m) they 
can be expressed in terms of modified Struve and Bessel functions. The cases to = 2 and m ~ 6 
are special in that remarkable simplifications occur: The asymptotic expansions of $2,5 and 
i>6,7 truncate at zeroth or first order in w"*, respectively, and these functions reduce to a simple 
exponential and a similar elementary function, respectively. This is what makes the calculation 
analytically feasible for to = 2, 6. 

To set up the RG for general to and d = d* (to) — e, one can introduce a renormalized field 
and renormalized coefficients a, r, p, and u via 

5p(y CFq^^"^ = nZpP , Uo CFq'^^'^ = jJi^ZuU , (11) 

'*More precisely, this applies to quantities retaining their physically significance for rri = 0, like r;L2(™=0) = f/ 
and vi^2{m=^) = ^- The series expansion coefficients of other exponents that are not required (nor have a known 
immediate physical meaning) in the m = theory, such as r)L4("i), may nevertheless have finite m -+ limits [6, 8]. 
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and determine the renormalization factors by minimal subtraction of poles. The residues of the 
Zs at two-loop order can be either analytically calculated or expressed in terms of single integrals 
whose integrands involve scaUng functions such as ^m,d* (m) • 

Let us illustrate this by means of the example of the two-loop graph which is 

oc G{x)^ in position space. To determine the Laurent expansion of this distribution, we apply 
it to a test function f{x\\ , x±), substituting Eq. (9) for G{x) and making the variable transfor- 
mation a;|| — !■ ■u = x\\/^/jr[. One is led to an integral of the form J^dx±^ x'^^~^ 9t.m{xi_) = 
(4e)~^5Q'^(0) -t- 0(e°). Upon evaluating this one sees that G{xY has poles oc A| 5{x)/e and 
a A_L 5{x) /e whose residues are proportional to 

POO POO 

/ dvv^'+^^m^'i^f and / dw""-^ ^m4'{^f , (12) 

respectively. Up to a simple factor given in Eq. (46) of Ref. [6], these are the integrals ja{fn) 
and j^im). The other integrals ja{fn) and Ju{fn) are of a similar nature but involve besides 
^rn,d* further scaling functions of the free theory like the one pertaining to the convolution 
(V||G * V||G)(cc) [= line with an insertion of (V||0)^]. 

The RG equations implied by the reparametrizations (11) and the /u-invariance of the bare 
theory can be exploited in a familiar manner to derive the scaling properties of correlation and 
vertex functions. For details the reader is referred to the original papers [5, 6], where the e 
expansions of the exponents are also utiUzed to obtain estimates for the numerical values of the 
critical exponents in the physical interesting cases with d = 3. The resulting values agree, in 
particular, remarkably well with the results of a very recent Monte Carlo investigation [34] of 
the ANNNI model, but also with some (though not all) experimental results and other theoretical 
estimates. 

So far we assumed perfect isotropy of the derivative terms of the Hamiltonian (3) in x\\- 
space. This is unrealistic in the case of lattice systems. Unless there is only one parallel direction 
(m = 1) or none, we should allow for a more general g| term, making the replacement 

m 

(A||<^)2 ^ WaT^^^i {didjcj>)dkdi4> = (A||</.)2 + w Y,{d!<Pf + ... (13) 

i=l 

in H, where T^"^; are the respective fourth-rank tensors compatible with the symmetry of the 
system. In the case of cubic synmietry, the only other term on the right-hand side besides the 
symmetric (first) contribution is the second one. To order e^, the dimensionless variable w; is 
found to be relevant at the isotropic w = fixed point (the coefficient of the term of its RG 
eigenexponentis positive) [50]. Unfortunately, the RG analysis for general w ^ that is required 
to decide whether a new fixed point exists to which the flow leads is rather compUcated and stiU 
remains to be completed. 



6 Concluding remarks 

Almost 25 years after the discovery of Lifshitz points and the introduction of the continuum mod- 
els (3), systematic field-theoretic RG analyses beyond one-loop order of the universality classes 
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these models represent have finally become feasible. This could open the way to accurate quanti- 
tative field-theory investigations and detailed comparisons with experimental results. Hopefully, 
this will also trigger further experimental work on the critical behavior at Lifshitz points. On 
the theoretical side, the progress reported here in applying the field-theoretic RG to the study of 
continuum models like (3) could lead to further insights regarding the question as to whether and 
under what conditions anisotropic scale invariance implies additional symmetries. 
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